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Abstract. We present an exact algebraic solution of a single graphene plane in 
transverse electric and perpendicular magnetic fields. The method presented gives 
both the eigen-values and the eigen-functions of the graphene plane. It is shown that 
the eigen-states of the problem can be casted in terms of coherent states, which appears 
& • in a natural way from the formalism. 
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■^j- ' 1. Introduction 

A major breakthrough in condensed matter physics took place when K. S. Novoselov et 
al. [1] , at Manchester University, UK, discovered an electric field effect in atomically 
thin carbon films. This electric field effect is characterized by the control of the 

O ■ 

electronic density in the films using a backgate setup. These atomically thin carbon 



X 



films were thought not to exist since long range order in two dimension can not occur. 
The system solves the apparent paradox by forming ripples. A single atomic layer 
of these thin carbon films is called graphene and its electric and magneto-electric 
properties triggered a new research field in condensed matter physics. The manufacture 
of graphene was followed by the production of other 2D crystals [2], which, however, 
have not the same exciting properties as graphene does. Applying high magnetic fields 
to a graphene sample, the Manchester group discovered that the quantization rule for 
the Hall conductivity is not the same one observes in the two-dimensional electron gas, 
being given instead by [3] 

a ™ = 4 ( n + I) J ' (1) 

with n an integer including zero. A confirmation of this result was independently 
obtained by Philip Kim's group [4J, at Colombia University, New York, USA. This new 
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Figure 1. Valence and conduction bands of electrons in graphene. The two bands 
touch each other in six points of the Brillouin zone, called Dirac points. 

quantum Hall effect was predicted by two groups working independently and using 
different methods [HE]. As explained by the two groups the new quantization rule for 
the Hall conductivity is a consequence of the dispersion relation of the electrons in the 
honeycomb lattice. This dispersion resembles the spectrum of ultra relativistic particles, 
i.e., the dispersion is that of particles of zero rest mass and an effective velocity of light. 
For graphene the effective velocity of light is Vp = c/300, with c the true velocity of 
light. 

For a qualitative description of the physics of graphene, both theoretical and 
experimental, see references [7] by Castro Neto et al, [8 J by M. I. Katsnelson, and 
[9] by Geim and Novoselov. 

In Figure [I] we show the energy dispersion of electrons in the honeycomb lattice. 
The spectrum shows a valence (lower) and a conduction (upper) bands. Since graphene 
has one electron per unit cell the valence band is completely filled and the properties 
of the system are determined by the nature of its spectrum close to the points where 
the valence and the conduction bands touch each other. These points are called Dirac 
points and are in number of six. In Figure [2] we show the spectrum close to the Dirac 
points. It is clear that the spectrum has conical shape of the form 

E = ±v F p, (2) 

where p is the magnitude of the momentum p around the Dirac point. 
A relativistic particle has an energy given by 

E = ^Jm 2 c A + p 2 c 2 , (3) 

and therefore an ultra-relativistic particle (m — > 0) has a spectrum given by 

E = cp. (4) 

It is clear from equation (J3J) that electrons in graphene, close to the Dirac points, 
have an energy dispersion with a formal equivalence to ultra-relativistic particles. As a 



Algebraic solution of a graphene layer in a transverse electric and perpendicular magnetic fields?) 




k... 



Figure 2. Valence and conduction bands close to one of the six Dirac points. It is 
clear that the bands can be approximated by a conical dispersion. 



consequence the quantum properties of the system has to be described by the massless 
(ultra-relativistic) Dirac equation in two plus one dimensions. 

We are interested in studying the spectrum of massless Dirac particles in the 
presence of a magnetic field perpendicular to the plane and an in plane homogeneous 
electric field, both static, a situation that occurs in the Hall effect. In the next section 
we present a full algebraic solution to this quantum problem. 

2. Algebraic solution 

2.1. Hamiltonian 

The massless Dirac equation in two plus one dimensions has the form 



where <7j, i = x,y,z, represents the Pauli spin matrices and pi, i = x,y, is the i 
component of the momentum operator in the position basis p = —ih'V. Since we 
are looking for stationary states we make the substitution ty(r,t) = ip(r)e~ tet ^ h . This 
substitution transforms the Dirac equation into an eigenvalue problem of the form 



The introduction of a magnetic field into a quantum mechanical problem is made by 
transforming the momentum operator according to the rule (minimal coupling) 



where A(r) is the vector potential and q is the charge of the particle. The magnetic 
field B is obtained from A using the relation B = V x A. There is a lot of freedom in 
choosing A and a common choice, known as Landau gauge, is A = (—By, 0, 0). Let us 
now assume that in addition to the magnetic field one has a homogeneous electric field, 




(5) 




(6) 




(7) 
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perpendicular to the magnetic field, and oriented along the y direction. This adds to 
the Hamiltonian a term of the form 

qVl = qSyl , (8) 

where V is the electric potential associated with the applied electric field E = (0, S,0) 
and 1 is the 2x2 unit matrix. 

Putting all together, the problem of a graphene layer in the presence of a magnetic 
field perpendicular to the layer and of an electric field parallel to the layer has its 
Hamiltonian, in the position basis, given by: 

H{x } y) = v F ( qSy/VF P*~iPy + <lBy\. (9) 
\ Px + iPy + qBy qtyj v F J 

The eigenproblem H(r)ip(r) = eip(r) can be further simplified by representing the 

eigenf unction ip(r) as 

V(r) = e ikx <j>{y) , (10) 

suggested by the translational invariance of Eq. (jUJ) along the x direction. Because we 
are dealing with electrons one has q = — e, with e > 0. Using equation (fTUI) in equation 
(JHJ) we obtain 

Next we perform a change of variables y = yls + l^k and d/dy = Isd/dy (corresponding 
to the introduction of the adimensional length y), with Is = \Jh/(eB) the magnetic 
length, and introduce the operators 

a =7f^ + 9 /^)' ( 12 ) 

a^ = ^(l-d/d§), (13) 

which satisfy the standard commutation relation [a, a'] — 1. Note that in Eqs. (jl2]) and 
( [TBI we have used the hat to distinguish between operators and their matrix elements 
in a given basis. 

Performing standard manipulations the Hamiltonian operator can be brought into 
a more transparent form 

~ = _ / eSl 2 B k + E B {a + at) E F a \ 

{ E F rf e£l 2 B k + E B (a + at) ) ' 1 ' 

with E F = sTlvphjln and Eb = eEls/V^- The eigenvalue equation one needs to solve 
has the form 

/ E B (a + a)) E F a \ ( |ai) \ _ / \a{) \ , . 

^ E F a^ E B {a + a^) ) \ \a 2 ) ) ~ € ° { \a 2 ) ) ' 1 

where eo = — (e + eSl B k). The eigenproblem is now in its most simplified form, with an 
effective Hamiltonian operator given by, 

n, _ ( E B (a + a^) E F a \ 

n ~{ E F a) E B (a + al) ^ 
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2.2. Diagonalization method 

Before going further with the diagonalization it is worth mentioning some properties of 
the Hamiltonian operator 7i, defined in Eq. ( fl6l) . First we define 7i as, 

fi = <rfa.= ( E *' + P " £ fA (17) 
\ —E F a J E B (a + a J ) J 

which, by definition, is an operator acting on the same Hilbert space as IK. Then we 
make the observation that both 'H + 'H and KK can be written as the 2x2 unit matrix 
times a simple operator (not a matrix operator) plus a 2 x 2 real matrix. The same 
holds for their linear combination, which enables us to write 

i2(n + n) + unn = ji + k, (is) 

for some simple operator J, some 2x2 real matrix K, and real /i and v. Now let the 
spinor be an eigenstate of K. Applying the left hand member of Eq. ( fl8i) to we 
obtain, 

\pt(H + H) + vHH\ |V) = {m + l#t + vcqH) |V) , (19) 
which means that if we chose // = eo and u = — 1 we reduce our problem to, 

(Jl + KM=e^>, (20) 

with 

j = (E 2 F - 2E 2 B )n - E 2 B (aa + o)a)) + 2e E B {a + a 1 ), (21) 
where n = a^a is the number operator, and 



K=\^-f E %°) (22, 



Equation fTSOj) is indeed simpler than our starting point, Eq. ffT5]) . We readily see 
that the spinor given by 



n I - (23) 
can be written as 

IV>> = 10) ( I ) , (24) 
where \<f>) is the eigenvector of the operator J and the spinor 

X R =( U ), (25) 



v 



is the right eigenvector of the eigenvalue problem 

K X R = X X R . (26) 

Nevertheless, there is one subtlety we must consider. Our simpler eigenproblem defined 
by Eq. (|2"U|) is such that their eigenvalues are eg, the squared eigenvalues of the original 
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problem given by Eq. (|T5|) . If e in Eq. (TT5|) and in Eq. (121711 have the same degeneracy 
it is guaranteed that eigenvectors of both problems are the same, and we can carry on 
the diagonalization with any of those equations. If, however, both +eo and — eo are 
eigenvalues of Ti, in Eq. (fT5l) . then our simple eigenproblem in Eq. (1201) will show an 
extra double degeneracy. This extra degeneracy must be handled carefully. Due to 
mixing the corresponding degenerate eigenvectors need not to be eigenstates of 7i, and 
these have to be found as particular linear combinations of the degenerated eigenvectors. 
It is instructive to switch off the electric field, for which Eq. ( {151) can be easily solved[5], 
and check whether the above problem shows up in the present case. 

2.3. Zero electric field case 

In the absence of the electric field we have £ = Eb = 0. As a consequence the operator J 
in Eq. ff2TT) is already in its diagonal form, being analogous to the ID harmonic oscillator: 
J\n) = Epn\n). In what regards the eigenproblem defined by Eq. (1261) for £ = 0, we 
can obtain the corresponding eigenvectors as 

*?=(J) and ^=(?)' (2T) 

with eigenvalues A + = E\ and A_ = 0, respectively. Recalling Eq. (1201) we get 
eg(n, ±) = E]?n + \±, where we recognize immediately the double degeneracy el(n, +) = 
6g(n + 1, — ) = Ep(n + 1). This degeneracy is in fact due to the presence of both +e 
and — eo as eigenvalues of H in the absence of electric field. Solving Eq. ( floT) directly 
for £ = gives [5] e = ±E F \/n + 1 in addition to the zero eigenvalue e = 0, with 
eigenvectors 

w± =( ± i!ti>) and w »=(io ))- < 28) 

respectively for nonzero and zero eigenvalues. Therefore we see that our method gives 
correctly |?/>)o, the only non-degenerate eigenvector, 

|V>>o=|0>X fl , (29) 

while \ip)± is given as the bonding and anti-bonding combination of degenerate 
eigenvectors an d \n + 

|V>)± = |n)x?±|n + l)xf ■ (30) 

As extra degeneracies due to the presence of a finite electric field are not to be expected, 
we will be able to identify any double degeneracy arising from Eq. (|20|) as a consequence 
of the presence of symmetrical eigenvalues ±eo in the original problem. 

2.4- Finite electric and magnetic fields 

With the above analysis in mind we proceed with the diagonalization of our problem 
in a finite electric and magnetic fields using Eq. (TSUI) . Let us start by solving Eq. f[2"S"j) . 
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The corresponding eigenvalues are given by, 

1 
2 

and as right eigenvectors we obtain 



X ± = -E 2 B + -[E 2 F ±E F ^E 2 -4E 2 B ) , (31) 



E 



B 



Ei 



(32) 



with 



1 



C ± = -(E P /\E B \±y/E? F /E? B -4) . (33) 



From equation fl3Tl) we see that E F and E B must satisfy the relation, 

E% > AE% , (34) 

if real eigenvalues are to be obtained. The meaning of this inequality is discussed later. 
Having solved the eigenproblem (1261) . the eigenproblem 

J\4>) = {el-X±)\4>) (35) 

remains to be solved. The solution of the eigenproblem (1350 is obtained in three steps. 
First we write the operator J as the sum of two terms, H\ + H 2 , given by 

E x = (E F - 2E B )n - E 2 B (aa + aV) , (36) 

and 

H 2 = 2e Q E B {a + a)). (37) 

As a second step we diagonalize the Hamiltonian H\ using the canonical transformation 

a) = cosh Urf — sinh t/7 , (38) 

and the corresponding Hermitian conjugated form for a. Replacing the a* and the a 
operators in ( 1361) one obtains 

H x = E 2 B 2 sinh U cosh U + (E F — 2E%) sinh 2 U + (7 t 7 t + 77) [-E|(cosh 2 U + sinh 2 [/) 
- (E 2 F - 2E 2 B ) sinh U cosh C/] + 7 t 7[(^| - 2E%) (cosh 2 C/ + sinh 2 U) 
+ 4E| sinh U cosh E7] . (39) 

In order for Hi to be diagonal we require that the coefficient multiplying the (7^7^ + 77) 
term should be null, leading to 

[-E% (cosh 2 U + sinh 2 U) - (E 2 F - 2E%) sinh U cosh U] = (40) 

which can be cast in the form 



^ 2U ^-e0Iei- < 41 > 

We note that since cosh U > for any value of C/ one must have sinh U < in order 
to satisfy equation (J4T]) . The result fT4Tl) together with cosh 2 [/ — sinh 2 f/ = 1 can be 
recasted in the form 

1 r 



sinh [/ = 

2 



l-( J E 2 -2i? 2 )/^ , (42) 
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and 

[l + (E 2 F - 2E 2 b )/uj] , (43) 



cosh 2 U = - 



leading to 



E 2 

sinht/cosht/ = - (44) 

uj 



with uj = J Ep — AEpE B and Ep > AE 2 B . Using the results for sinhf/ and cosh ?7 one 
can write the piece Hi of the full Hamiltonian as 

Hi = ±[u>- (E 2 F - 2E 2 B )} + uy^ = d + w T t 7 , (45) 

with Ci = [uj — (Ep — 2E B )]/2. The piece H 2 has now the form 

H 2 = 2e E B (coshU- sinh f/)( 7 t + 7 ) = C 2 (^ + j) , (46) 

where C 2 = 2e E B (coshU — sinh £7). The third step requires the diagonalization of 
Hi + H 2 in the new form, written in terms of the 7 — operators; this is accomplished by 
the transformation 7 t = (ft + Z (with Z = —C 2 /uj), leading to 

Hi + H 2 = Ci-^- + (47) 

UJ 

which has the desired diagonalized form. The eigenenergies of the Hamiltonian (14711 
have the form, 

uJ n = Ci-^- + ujn = huj-(E 2 p- 277/1)] - + uju , (48) 

uj Z uj 

and the ground state obeys the relation 

/3|0;/?)=0<^ 7 |0;/3) = Z|0;/3>. (49) 

One should note that the state |0; (3) is an eigenstate the 7 operator with eigenvalue Z\ 
it is therefore said that |0; (3) is a coherent state of the operator 7 . This last result allows 
us to write the vacuum of the (3 operators in terms of the vacuum of the 7 operators as 

|0;/3) = e^ t |0; 7 ), (50) 

and any eigenstate is written in terms of the /3-operators as 

|n; ft = 4=^)10; P) = 4^ (7 f - Z) n |0; 7 ) . (51) 



The eigenenergies 6q of our simpler eigenproblem defined in Eq. (1201) are obtained from 
[see Eq. ( 1551) ] 6q — A± = uj n , leading to 

e 2 (n,±) = ^[n + (l±l)/2]. (52) 
E F 

The double degeneracy €q(u, +) = ej^n+l, — ) for n ^ is immediately recognized. From 
the analysis we have made in Sees. 12.21 and 12.31 it is now obvious that this degeneracy 
signals the presence of both solutions ±eo in the original problem [Eq. (fT5l) ]. Moreover, 
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as the eigenvectors of the finite electric field problem have to equal Eq. 
we arrive at the following solution, 



e(n, ±) = —e£l 2 B k =F 



E 



1/2 



-Vn + 1 



when £ — > 



(53) 



with eigenvectors given by 



E 



B 



E f 



-\n;P)^T\n + l;P)y/UZ 
|n;/3)/V^±|n + l;/3)/V^I 



(54) 



where C± is defined in Eq. 



In addition there is a single non-degenerate solution 
0, which gives e = —e£l B k, and has as eigenvector 



E 



B 



E f 



|0;/3). 



(55) 



This concludes our solution. The eigen-values (|53|) agree with those obtained by 
Lukose et al. [ID]. These authors solved the problem directly in the position basis by 
transforming the original problem, by means of a Lorentz boost transformation, into a 
case where the electric field is null. 



2.5. Physical interpretation 

The standard 2D electron gas pierced by a perpendicular magnetic field is known after 
Landau [UJ to have a spectrum given by, 

e(n)=hujc(n + -), n = 0,l,2,... (56) 

in complete analogy with the quantum harmonic oscillator, where Uq = \eB\/m is the 
cyclotron frequency for electrons with mass m. The so-called Landau levels are equally 
spaced with level separation hu> c , which increases linearly with B. An in-plane electric 
field is easily handled by the transformation — + e£lB/(huj c \/2), whose major 
consequence is a shift of the entire spectrum, 

e(n) = -e£l 2 B k-^^ + htu c (n + ^), n = 0, 1, 2, . . . (57) 

with no change for the cyclotron frequency. 

Landau levels in graphene are completely different from Landau levels in standard 
2D electron gas. As mentioned in Sec. 12.31 graphene's spectrum in perpendicular 
magnetic field is given by [S], 

e{n) = ±Tiu cy fn, n = 0,1,2,... (58) 

where the cyclotron frequency is uj c = \/2vf/Ib- Two major differences become apparent 
when comparing Eqs. ( 1581) and ( |56|) . Firstly, Landau level spacing in graphene is not 
constant due to the square root in Eq. (1581 . Secondly, the standard 2D electron gas 
has ujq oc B whereas graphene shows Cue oc \fB. These dissimilarities are a direct 
consequence of the effective ultra relativistic nature of the quasi-particles in graphene. 
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Figure 3. First ten Landau levels as function of The quantity E c is the 

cyclotron energy huj c . The momentum k was chosen to be zero. 

As regards the presence of an in-plane electric field and perpendicular magnetic 
field in graphene we have shown that Landau levels are given by Eq. (153]) . which can be 
cast in the form, 

e(n) = -eS l%k T- hn cv ^n, n= 0,1,2,... (59) 
where the new cyclotron frequency reads 

n c = V2 V -^[l-£ 2 /(B 2 vl)f. (60) 

Thus, unlike the usual 2D electron gas, graphene's cyclotron frequency is renormalized 
by the electric field, as can be seen in Eq. (f60T) . which, of course, reduces to lj c in the 
limit £ — > 0. More important though is the fact that \S\ < f^l-Bl must be realized if 
Q c is to be real. This last inequality is exactly the same expressed in Eq. and its 
meaning is now unveiled. As £ approaches vpB from below Qc becomes smaller and 
smaller, and Landau levels become closer and closer. Eventually, the electric field is 
such that \£\ = vp\B\, which implies Qc — 0, and consequent collapse of Landau levels. 
In Fig. [3] we show the first ten Landau levels (for positive and negative energies) as 
function of \£ \/vp\B\\ the collapse of the Landau levels is clear. For \£\ > Vp\B\ the 
present solution is not valid. 

3. Concluding remarks 

The problem of a single graphene plane in transverse electric and perpendicular magnetic 
fields assembles in a simple way several algebraic methods of diagonalizing bilinear 
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problems. The matrix form of the Dirac Hamiltonian - the low energy effective 
Hamiltonian for graphene - calls for several non-standard manipulations where canonical 
transformations and the concept of coherent state appear in a natural way. Furthermore, 
alike the standard 2D electron gas pierced by a magnetic field, an additional in-plane 
electric field in graphene induces cyclotron frequency renormalization. Moreover, when 
the electric field equals the critical value v F B Landau level collapse is observed. 
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